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Abstract 

A  Fortran  computer  program  LOWTHRM  is  described  for 
calculating  nuclear  thermal  fluence  incident  upon  a  target 
area.  Atmospheric  transmissivity  factors  in  the  spectral 
region  0.25  to  28.5  microns  are  determined  through  use  of 
the  LOWTRAN 5  computer  code.  The  program  provides  a  choice 
of  six  model  atmospheres  covering  seasonal  and  latitudinal 
variations  from  sea  level  to  100  km,  eight  haze  models,  and 
accounts  for  molecular  absorption,  molecular  scattering,  and 
aerosol  extinction.  Atmospheric  refraction,  earth  curvature 
effects,  thermal  scattering,  and  thermal  ground  reflection 
contributions  are  included. 
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LOWTHRM: 

A  THERMAL  FLUENCE  CODE 

I.  Introduction 

Background 

One  of  the  unique  characteristics  of  a  nuclear  detona¬ 
tion  is  the  release  of  an  enormous  amount  of  energy  in  a  small 
mass  of  material  which  ultimately  results  in  the  formation  of 
a  fireball.  This  fireball  reradiates  part  of  its  energy  as 
thermal  radiation  distributed  over  the  ultraviolet,  visible, 
and  infrared  regions  of  the  spectrum.  It  is  this  thermal 
radiation  that  can  cause  a  rapid  heat  buildup  producing  burns, 
charring,  and  possibly  ignition  of  exposed  materials  (Ref.  9). 

The  energy  flux  of  thermal  radiation  from  a  nuclear 
detonation  decreases  rapidly  with  increasing  distance  from  the 
burst  point.  This  is  largely  due  to  spherical  divergence,  but 
atmospheric  factors  such  as  absorption  and  scattering  are  also 
important.  In  fact,  these  atmospheric  factors  change  the  radi¬ 
ation  incident  upon  a  target  from  a  direct  fluence  to  more  of 
a  diffuse  one.  Therefore,  a  target  receives  radiation  that  has 
been  scattered  toward  it  as  well  as  the  radiation  that  has 
been  directly  transmitted  from  burst  to  receiver.  Since  ther¬ 
mal  radiation  calculations  are  important  in  predicting  the 
survivability/vulnerability  of  Air  Force  systems,  a  model  for 
calculating  thermal  fluence  is  essential  to  Air  Force  planning. 
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Purpose 

The  purpose  of  this  thesis  is  to  develop  a  thermal  fluence 
computer  code.  It  will  determine  contributions  from  direct, 
scattered,  and  ground  reflected  sources  incident  upon  a  target 
area.  Values  for  "exact”  atmospheric  transmission  of  radia¬ 
tion  will  be  based  upon  the  L0WTRAN5  Fortran  computer  code 
currently  in  use  by  the  Air  Force.  To  determine  the  validity 
of  the  developed  code,  results  obtained  will  be  compared  against 
several  models  available  for  similar  computations. 

Scope 

In  order  to  alter  the  L0WTRAN5  code  to  produce  thermal 
fluences,  a  thorough  study  of  the  code  and  its  variables  must 
first  be  accomplished.  Next,  routines  to  calculate  the  weapon 
thermal  efficiency  and  to  weight  the  source  spectrum  must  be 
added.  Numerical  integration  techniques  are  then  used  to 
determine  the  direct,  scattered,  and  ground  reflected  fluences 
incident  on  a  receiver  at  various  orientations.  Finally,  a 
User's  Guide  is  included  to  document  options  available  from 
the  completed  code. 

Development  of  the  Report 

In  Chapter  II,  a  short  overview  of  the  theory  of  thermal 
radiation  from  an  air  burst  is  presented.  In  this  chapter 
the  source  is  presented  with  a  transmission  model  based  on 
diffusion  theory.  The  energy  range  of  the  thermal  photons 
and  mechanisms  of  interaction  with  the  atmosphere  are  examined. 
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Chapter  III  discusses  the  L0WTRAN5  code  upon  which  the  avail¬ 
able  atmospheric  models  are  based.  In  Chapter  IV  the  LOWTHRM 
computer  model  is  presented.  It  includes  discussion  of  the 
thermal  pulse  and  the  unreacted,  scattered,  and  ground  reflected 
fluences.  Chapter  V  compares  results  obtained  from  the  fin¬ 
ished  code  with  other  available  models.  Appendix  A  is  the 
LOWTHRM  User's  Guide  which  may  be  reproduced  and  used  as  an 
aid  in  running  the  code.  Appendix  B  contains  sample  output 
for  several  example  cases. 


II.  Theory  of  Thermal  Radiation  Phenomenon 

Source 

The  source  of  thermal  radiation  in  a  nuclear  detonation 
is  the  fireball.  For  an  air  burst,  the  thermal  radiation  is 
emitted  in  two  pulses  which  correspond  to  the  apparent  rise 
and  fall  of  the  fireball  surface  temperature.  The  first  pulse 
is  of  very  short  duration  and  contains  only  about  one  percent 
of  the  total  thermal  radiation  emitted.  The  second  pulse  con¬ 
tains  the  bulk  of  the  thermal  radiation  and  may  last  for  sever¬ 
al  seconds  depending  upon  the  weapon  yield  (Ref.  9).  The  two 
pulses  are  illustrated  in  Figure  1. 


TIME  AFTER  EXPLOSION 
(RELATIVE  SCALE) 


Fig.  1.  Emission  of  Thermal  Radiation  in  Two  Pulses  in 
an  Air  Burst  (Ref.  9:41) 
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Because  of  its  low  yield,  the  first  pulse  can  be  neglected  in 
thermal  radiation  calculations. 

In  order  to  construct  an  approximate  source  model,  some 
simplifying  assumptions  must  be  made.  The  nuclear  burst  is 
considered  to  be  a  point  isotropic  source  localized  at  the 
point  of  detonation.  It  is  also  assumed  to  emit  radiation 
according  to  a  Planckian  blackbody  distribution.  Although 
studies  made  in  various  weapons  tests  have  shown  that  the  fire¬ 
ball  does  not  behave  exactly  like  a  blackbody  radiator,  this 
distribution  is  close  to  the  actual  weapon  spectrum  (Ref.  9:305). 
To  the  extent  of  this  report,  the  deviation  from  an  ideal 
Planckian  spectrum  does  not  contribute  a  serious  error  to  the 
calculations  of  thermal  fluence. 

Transmission 

The  diffusion  equation  is  used  to  develop  the  theory  of 
transmission  of  thermal  radiation.  Assuming  a  homogenous 
atmosphere  and  all  scatter  as  removal,  the  thermal  radiation 
can  be  assumed  to  act  like  streaming  particles.  Taking  into 
account  the  point  source,  the  fluence,  F(r),  is  a  function 
of  the  distance  from  the  source.  The  diffusion  equation  is 

d^r-)  +  |  F(r)  +  yF(r)  =  0  (1) 

where  r  is  the  distance  from  the  source  and  y  is  the  total 
cross  section.  Rearranging  Equation  (1)  and  integrating 
from  the  fireball  radius,  r^r  to  the  receiver,  R,  gives 
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which  is 


,R  dF  (r) 

rfb  F(r) 


2  ,  fR 

jdr  -  v/ 
r  rfb 


dr 


F(R) 


(r2fb  F(rfb>  )exp  (-y(R-rfb)) 
R2 


(2) 


(3) 


Applying  boundary  conditions 

lim  (4irr2F  (r) )  =  Y  (4) 

r->rfb 

where  Y  is  the  thermal  yield 

F(rfb)  =  - 7 —  (5) 

4irr2fb 

Since  the  point  source  assumes  rfb  is  zero,  substituting 
Equation  (5)  into  Equation  (3)  yields 

f (r)  *  <.n«q?<aw.  (6) 

4irR2 

The  expression  exp(-yR)  is  the  atmospheric  transmittance,  T. 
This  transmittance  is  a  function  of  many  variables.  However, 
if  the  effects  of  the  atmosphere  are  ignored.  Equation  (6) 
reduces  to 

F  =  — 2L.  (7) 

4ttR2 

where  F  is  the  thermal  fluence  (energy  per  unit  area) ,  Y  is 
the  thermal  yield  (energy) ,  and  R  is  the  path  length  between 
the  burst  and  receiver  (distance) .  The  only  attenuation  in 
this  case  is  spherical  divergence  given  by  4ttR2  (Ref.  12:4-6). 
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Spectrum  from  an  Air  Burst 

The  range  of  thermal  photon  energies  from  the  source 
is  assumed  to  correspond  to  a  6000° K  Planckian  distribution 
(Ref.  9:307).  Using  Planckian  tables  (Ref. 23),  the  range  from 
0.25  to  28.5  microns  contains  approximately  99  percent  of  the 
thermal  radiation  emitted.  This  range  can  be  divided  into 
three  broad  bands;  ultraviolet,  visible,  and  infrared.  Table 
I  gives  the  band  limits  in  terms  of  wavelength  and  frequency. 

Table  I 

Broad  Band  Structure  of  Thermal 
Spectrum  in  an  Air  Burst 


BAND  WAVELENGTH  FREQUENCY 

(microns)  (Hertz/Sec) 


Ultraviolet 

0.25  - 

0.4 

1.2 

X 

1015  -  7.5 

X 

1014 

Visible 

0.4  - 

0.7 

7.5 

X 

1014  -  4.3 

X 

1014 

Infrared 

0.7  - 

28.5 

4.3 

X 

1014  -  1.1 

X 

10  1 3 

Photon  Interaction  with  the  Atmosphere 

Aerosol .  Of  course  the  transmission  effects  of  the 
atmosphere  cannot  be  ignored  since  this  would  predict  fluences 
that  are  higher  than  can  realistically  be  expected.  In  the 
visible  and  infrared  portions  of  the  spectrum,  scattering 
effects  drastically  change  atmospheric  transmission.  If  the 
atmosphere  were  composed  solely  of  air  molecules,  the  know¬ 
ledge  of  their  scattering  laws  would  be  sufficient  for 
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transmission  calculations.  However,  the  atmosphere  also  in¬ 
cludes  haze  aerosol  which  is  subject  to  other  scattering  laws. 
This  aerosol  is  composed  of  all  solid  or  liquid  particles 
suspended  in  the  atmosphere  with  sizes  ranging  from  about  0.03 
to  100  microns.  Although  aerosol  concentrations  are  orders 
of  magnitude  less  than  molecular  concentrations,  the  aerosols 
are  very  important  in  atmospheric  transmission  of  radiation. 

These  particles  scatter  and  absorb  radiation  and  affect  the 
processes  of  atmospheric  condensation  resulting  in  the  forma¬ 
tion  of  clouds,  snow,  rain,  and  fog  (Ref.  2).  This  scattering 
of  thermal  radiation  is  generally  treated  according  to  the 
classical  theory  developed  by  Gustav  Mie  (Ref.  22) . 

Molecular  Scatter.  At  altitudes  greater  than  five  kilo¬ 
meters  the  aerosol  concentration  approaches  zero  and  scattering 
is  principally  due  to  the  molecular  constituents  of  the  atmos¬ 
phere.  Scattering  from  such  particles  is  termed  Rayleigh 
scattering  after  Lord  Rayleigh,  who  in  1871  developed  the 
initial  theory  for  scattering  of  radiation  by  air  molecules. 
Rayleigh's  law  is  based  upon  the  assumption  that  the  air  mole¬ 
cules  are  homogenous  spherical  particles  of  equal  size  whose 
radii  are  much  smaller  than  the  wavelength  of  the  incident 
electromagnetic  radiation.  This  is  clearly  the  case  since 
atmospheric  molecules  have  dimensions  of  10” 6  to  10” 8  microns, 
much  smaller  than  the  0.25  to  28.5  micron  range  under  study. 

After  these  particles  are  hit  by  the  incident  radiation,  they 
become  a  source  of  electromagnetic  radiation  of  a  specific 

wavelength  transmitted  in  all  directions  whose  field  distribution 
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is  described  by  Maxwell's  equations.  Rayleigh's  theoretical 
approach  only  covers  primary  scattering  and  reveals  that  the 
scattering  coefficient  due  to  air  molecules  varies  inversely 
with  the  fourth  power  of  the  wavelength  (Ref.  5) . 

Absorption 

Radiation  is  scattered  by  all  of  the  atmospheric  gaseous 
molecules,  but  is  selectively  absorbed.  Minor  constituents 
of  the  atmosphere  such  as  water  vapor,  carbon  dioxide,  ozone, 
nitrous  oxide,  carbon  monoxide,  and  methane  are  dominant 
factors  in  the  absorption  process.  The  two  major  atmospheric 
constituents,  nitrogen  and  oxygen,  are  minor  absorbers  when 
compared  to  the  preceeding  list.  Atmospheric  ozone  is  par¬ 
ticularly  effective  in  absorbing  ultraviolet  radiation  while 
the  infrared  is  susceptible  to  water  vapor  and  carbon  dioxide. 
For  a  detailed  explanation  of  atmospheric  absorption  see 


Reference  4 . 
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III.  L0WTRAN5  Computer  Code 

In  order  to  calculate  thermal  radiation  incident  on  a 
target  area,  the  transmissivity  of  the  atmosphere  must  first 
be  accurately  determined.  The  Fortran  computer  code,  L0WTRAN5, 
is  designed  to  make  this  calculation  for  a  given  atmospheric 
path  over  the  frequency  range  of  0.25  to  28.5  microns  on  a 
linear  wavenumber  scale.  A  choice  of  six  model  atmospheres 
is  given  with  an  option  for  a  seventh  model  which  can  be  in¬ 
serted  as  a  set  of  radiosonde  data.  Aerosol  attenuation  is 
calculated  for  a  given  visual  range  based  on  an  interpolation/ 
extrapolation  scheme  using  any  two  of  eight  aerosol  models 
provided. 

For  horizontal  path  transmittance  calculations  under 
nonstandard  conditions,  the  user  can  specify  meteorological 
data  for  use  in  the  program. 

Theory 

In  determining  the  transmittance  of  the  atmosphere, 
LOWTRAN5  follows  almost  exactly  the  procedures  outlined  by 
McClatchey  et  al  (Ref. 13).  The  main  assumptions  made  are 
that  the  atmosphere  can  be  represented  by  a  33-layer  model, 
and  that  the  average  transmittance  over  a  20  cm-1  interval 
can  be  represented  by  a  single  parameter  model  of  the  form 

T  =  f  (CW)  .  (8) 

C  is  a  wavelength  dependent  absorption  coefficient  and 
W  is  an  "equivalent  absorber  amount"  for  the  atmospheric  path. 
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This  atmospheric  path  is  defined  in  terms  of  the  pressure, 
p(z),  temperature,  T(z),  concentration  of  absorber,  L,  and 
an  empirical  constant  as  follows: 


W  =  L< 


P(z) 

P 


/  T  . 
/  T (z)  ' 


n 


(9) 


If  Equation  (9)  is  substituted  into  Equation  (8)  and  n 
is  set  equal  to  zero  and  unity,  respectively.  Equation  (8) 
reverts  to  the  well  known  weak  line  and  strong  line  approxi¬ 
mations  common  to  most  band  models  (Ref. 18). 

The  form  of  the  function  f  and  parameter  n  was  determined 
empirically  using  both  laboratory  transmittance  data  and 
available  molecular  line  constants.  In  both  cases,  the  trans¬ 
mittance  was  degraded  in  resolution  to  20  cm-1  throughout  the 
entire  spectral  range. 


Model  Atmospheres 

The  altitude,  pressure,  temperature,  water  vapor  density, 
and  ozone  density  for  the  U.S.  Standard  Atmosphere  and  five 
seasonal  model  atmospheres  are  provided  as  basic  input  data 
to  LOWTRAN5.  Also  provided  are  the  number  of  aerosol  parti¬ 
cles  per  cubic  centimeter  for  eight  haze  models .  The  model 
atmospheres  correspond  to  the  1962  U.S.  Standard  Atmosphere, 
Tropical  (15°N)  Atmosphere,  Midlatitude  Summer  (45°N,  July) , 
Midlatitude  Winter  (45°N,  January),  Subarctic  Summer  (60°N, 
July)  and  Subarctic  Winter  (60°N,  January).  The  different 
models  are  digitized  in  1  kilometer  steps  from  0  to  25  km,  5 
km  steps  from  25  to  50  km,  then  at  70  km  and  100  km  directly 

as  given  by  McClatchey  (Ref. 17). 
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The  water  vapor  and  ozone  altitude  profiles  added  to  the 
1962  U.S.  Standard  Atmosphere  by  McClatchey  were  obtained 
from  Sissenwine  (Ref. 21)  and  Herring  (Ref. 10)  respectively, 
and  correspond  to  mean  annual  values.  The  water  vapor  densi¬ 
ties  for  the  1962  U.S.  Standard  Atmosphere  correspond  to 
relative  humidities  of  approximately  50  percent  for  altitudes 
up  to  10  kilometers.  However,  the  relative  humidity  values 
for  the  supplementary  models  tend  to  decrease  with  altitude 
from  approximately  80  percent  at  sea  level  to  30  percent  at 
10  kilometer  altitude. 

In  addition  to  the  model  atmospheres  provided,  the  user 
has  the  option  of  inserting  his  own  model  atmosphere,  or  of 
building  another  model  by  combining  various  parts  of  the  six 
standard  models  (Ref. 18). 

Atmospheric  Constituents 

Atmospheric  Gases .  LOWTRAN5  assumes  that  mixing  ratios 
of  the  gases,  carbon  dioxide,  nitrous  oxide,  methane,  carbon 
monoxide,  nitrogen,  and  oxygen  remain  constant  at  all  alti¬ 
tudes  at  the  values  listed  in  Table  II.  These  gases  as  a 
whole,  with  the  exception  of  nitrogen,  are  referred  to  as 
the  uniformly  mixed  gases. 

Absorption  coefficients  for  water  vapor,  ozone,  and  the 
combined  effects  of  the  uniformly  mixed  gases  were  digital¬ 
ized  from  spectral  curves  by  McClatchey  and  are  included  as 
data  for  LOWTRAN  5  (Ref. 13). 
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Table  II 


Mixing  Ratios  for  Atmospheric 
Gases  in  Parts  per  Million 


CARBON  DIOXIDE 

330 

NITROUS  OXIDE 

0.28 

METHANE 

1.6 

CARBON  MONOXIDE 

0.075 

NITROGEN 

7.905x10s 

OXYGEN 

2.095x10s 

Continuum  Absorption.  Absorption  coefficients  for  the 
water  vapor  continuum  are  based  on  measurements  of  Burch 
(Ref.  3),  McCoy  and  Rensch  (Ref. 14),  and  Bignell  (Ref.  1). 

The  continuum  due  to  collision  induced  absorption  by 
nitrogen  is  included  based  on  the  measurements  of  Reddy  and 
Cho  (Ref. 16)  and  Shapiro  and  Gush  (Ref. 20). 

In  all  cases  the  transmittance  due  to  continuum  absorp¬ 
tion  is  assumed  to  follow  a  simple  attenuation  law. 

Molecular  Scattering.  The  absorption  coefficient  due 
to  molecular  scattering  in  introduced  into  LOWTRAN5  via  the 
following  expression: 

M.S.  =  (9 . 807x10” z  0 )  (W4-0117)  km-1  (10) 

where  W  is  in  wavenumbers  (cm-1). 

The  above  expression  was  obtained  as  a  best  fit  to 
molecular  scattering  coefficients  published  by  Penndorf  (Ref. 15). 
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Aerosol  Models.  Eight  aerosol  models  are  incorporated 
into  LOWTRAN 5 .  These  include  Rural  (23km  visibility).  Rural 
(5km  visibility).  Maritime  (23km),  Maritime  (5km),  Urban  (5km), 
Tropospheric  (50km),  Fogl  (0.2km),  and  Fog2  (0.5km).  Although 
specific  visibilities  are  listed  for  each  model,  an  aerosol 
attenuation  for  any  visual  range  can  be  calculated  using  an 
interpolation/extrapolation  procedure . 

The  Maritime,  Urban,  Rural,  and  Average  Continental 
Aerosol  Models  are  all  boundary  layer  models.  That  is,  they 
apply  to  the  first  few  kilometers  of  the  atmosphere.  The 
Tropospheric  Model  was  developed  primarily  for  use  in  the 
troposphere  above  the  boundary  layer.  However,  it  can  be 
used  for  calculations  near  ground  level  for  particularly 
clear  and  calm  conditions  in  pollution  free  areas  with  visi¬ 
bilities  greater  than  thirty  to  forty  kilometers. 

The  aerosol  models  are  based  on  the  following  assumptions. 
A  particle  size  distribution  similar  to  Deirmendjian 1 s  Haze 
Model  C,  but  where  the  large  particle  radius  cutoff  has  been 
extended  to  one  hundred  microns  (Ref.  6) .  The  particle  size 
distribution  is  assumed  to  remain  constant  with  altitude 
(Ref. 18).  The  variation  of  aerosol  number  density  with  alti¬ 
tude  is  assumed  to  be  the  same  as  given  by  McClatchey  for  the 
twenty-three  kilometer  visual  range  model.  The  latter  aerosol 
number  densities  were  adjusted  to  give  extinction  coefficients 
at  a  wavelength  of  0.55  microns  that  correspond  to  those 
obtained  by  Elterman  at  each  altitude  (Ref.  7) .  The  variation 
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of  aerosol  refractive  index  with  wavelength  has  been  obtained 
from  measurements  by  Volz,  who  found  that  aerosols  are  com¬ 
posed  of  water-soluble  substances  as  well  as  dust-like  material 
(Ref. 25) . 

Aerosol  extinction  and  aerosol  absorption  values  were 
calculated  based  on  single  scattering  Mie  theory  using  the 
above  aerosol  size  distribution  and  refractive  index  values. 

The  Tropospheric  Aerosol  Model  is  the  same  as  the  small 
particle  portion  of  the  Rural  and  Urban  models.  The  larger 
particles  are  lost  at  a  higher  rate  than  the  small  ones,  and 
above  the  boundary  layer  they  are  not  replaced  by  turbulent 
mixing  from  the  surface.  The  continental  component  of  the 
Maritime  Model  also  is  the  same  as  the  small  particle  portion 
of  the  Rural  Model  for  analgous  reasons. 

The  aerosol  models  also  differ  in  their  composition  and 
the  corresponding  variation  of  refractive  index  with  wave¬ 
length.  The  Rural  Model  is  assumed  to  be  a  mixture  of  seventy 
percent  water-soluble  aerosols  and  thirty  percent  dust-like 
aerosol . 

The  Maritime  Model  is  composed  of  a  mixture  of  aerosols 
of  oceanic  and  continental  origins.  The  oceanic  aerosols 
are  produced  primarily  by  the  sea  spray  and  are  assumed  to 
be  a  solution  of  sea  salts  in  water.  The  continental  com¬ 
ponent  has  the  same  composition  as  the  Rural  Model.  While 
the  proportions  and  nature  of  the  two  components  of  the  Mari¬ 
time  aerosol  will  vary  geographically,  there  is  insufficient 


data  to  meaningfully  model  these  variations.  For  simplicity, 
the  oceanic  component  is  taken  as  contributing  seventy-five 
percent  of  the  extinction  at  0.55  microns,  which  yields  a 
model  which  is  consistent  with  measurements  in  a  number  of 
different  maritime  locations  (Ref. 19). 

The  Urban  Model  is  similar  to  the  Rural  Model,  but  with 
an  addition  of  soot-like  aerosols  such  that  the  mixture  is 
thirty-five  percent  soot-like  aerosols  and  sixty-five  percent 
rural  aerosols.  The  Tropospheric  Model  is  assumed  to  have 
the  same  composition  as  the  rural  aerosols  (Ref. 19). 

The  characteristics  of  the  different  aerosol  models,  for 
the  lower  atmosphere,  are  summarized  in  Table  III.  The  size 
distributions  are  represented  by  one  or  the  sum  of  two  log¬ 
normal  distributions: 

dN(r)  2  Nj  (log  r  -  log  r^) 2 

-  =  I( - )  exp  ( - )  (11) 

dr  i=l  £n  (10)  ra^/TiF  2ajL2 

The  choices  of  N0  in  Table  II  correspond  to  one  particle  per 
cubic  centimeter.  The  actual  size  distributions  are  renor¬ 
malized  to  give  the  correct  extinction  coefficients  and  models 
being  used  (Ref. 19). 

Aerosol  Interpolation/Extrapolation  Scheme.  The  total 
extinction  coefficient,  s (T) ,  at  0.55  microns  is  inversely 
proportional  to  visual  range,  VIS,  and  can  be  written  as 
follows : 

3.91 

s(T)  =  s(a)  +  s (m)  =  -  (12) 

VIS 
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Table  III 

Characteristics  of  the  Aerosol 
Models  of  the  Lower  Atmosphere 


AEROSOL  MODEL 

SIZE 

DISTRIBUTION 

TYPE 

i 

Ni 

ri 

ai 

Rural 

1 

0.9999975 

0.005  ym 

0.475 

Water-Solubles  and 

2 

0.0000025 

0.5  ym 

0.475 

Dust-Like 

Urban 

1 

0.9999975 

0.005  ym 

0.475 

Rural  Aerosol 

Mixture  and 

2 

0.0000025 

0.5  ym 

0.475 

Soot-Like 

Maritime 

Continental 

Origin 

1.0 

0.005  ym 

0.475 

Rural  Aerosol 

Mixture 

Marine 

Origin 

1.0 

0.3  ym 

0.4 

Sea  Salt  Solution 
in  Water 

Tropospheric 

1.0 

0.005  ym 

0.475 

Rural  Aerosol 

Mixture 

This  assumes  a  two  percent  contrast  threshold  where  the  "a" 
and  "m"  refer  to  the  aerosol  and  molecular  components  respec¬ 
tively.  The  aerosol  extinction  coefficient  can  then  be 
written  as: 

s (a)  =  -  s (m)  (13) 

Since  the  aerosol  extinction  coefficient  s (a)  is  directly 
proportional  to  the  aerosol  number  density  N(z),  it  can  be 
written 

N(Z)  =  +  b ( Z )  (14) 

where  a(Z)  and  b(Z)  are  constants  for  a  given  altitude  Z. 

It  will  be  noted  that  b(Z)  is  proportional  to  the  molecular 
scattering  coefficient  at  0.55  microns  at  altitude  Z,  where 
a  molecular  absorption  has  been  assumed  negligible  at  \  =  0.55. 

The  above  equation  forms  the  basis  for  the  interpolation/ 
extrapolation  procedure  used  in  LOWTRAN 5  to  determine  the 
aerosol  attenuation  at  any  given  visual  range. 

The  coefficients  a  and  b  are  determined  from  the  above 
equation  at  each  altitude  using  the  two  aerosol  models,  that 
is, 

a(Z)  =  (N5  (Z)  -  N23(Z)>  /  d/5  -  1/23)  (15) 

b(Z)  =  (N23(Z)/5  -  N5(Z)/23)  /  (1/5  -  1/23)  (16) 
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where  N5  and  N23  refer  to  the  number  densities  for  the  five 
kilometer  and  twenty-three  kilometer  visual  ranges.  The  above 
procedure  is  used  only  in  the  lower  five  kilometers  of  the 
atmosphere  since  the  two  aerosol  models  are  identical  above 
five  kilometer  altitude  (Ref. 19). 
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IV.  Computer  Model  for  Thermal 
Fluence  Calculations 

Because  of  the  difficulties  involved  in  what  is  essen¬ 
tially  a  three-dimensional  transport  problem,  certain  simpli¬ 
fying  assumptions  had  to  be  made  in  order  to  produce  a  fast¬ 
running  thermal  fluence  code.  LOWTHRM,  in  its  final  form, 
is  the  result  of  several  preceeding  "generations"  of  the 
original  computer  model.  As  different  ideas  were  incorporated 
or  discarded,  the  code  changed  from  basically  an  atmospheric 
transmittance  model  to  an  efficient  thermal  fluence  code. 

This  chapter  introduces  the  thermal  source  and  geometrical 
models  used  in  the  program  along  with  the  development  of  the 
normalized  scatter  integral. 

Thermal  Pulse 

As  stated  earlier  in  this  report,  the  nuclear  detonation 
is  assumed  to  be  an  isotropic,  spherical  point  source  at  the 
burst  location.  Since  fluence  calculations  will  be  for  dis¬ 
tances  relatively  far  from  the  burst  point,  this  approximation 
should  yield  reasonable  results  (Ref.  24).  A  6000°K  blackbody 
is  generally  used  to  approximate  a  nuclear  airburst.  For  low 
altitudes  this  blackbody  is  not  a  bad  approximation.  As  the 
burst  altitude  is  increased,  however,  the  spectrum  shifts 
toward  the  ultraviolet  wavelengths.  To  compensate  for  this 
shift,  an  interpolation  routine  was  developed  by  Johnson  (Ref. 
2)  to  raise  the  blackbody  temperature  as  burst  height  is 
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increased.  The  6000°K  blackhody  is  used  for  all  burst  heights 
15  km  or  less.  A  7000°K  blackbody  is  the  model  for  a  burst 
at  25  km.  For  bursts  between  15  and  25  km,  linear  interpola¬ 
tion  is  used  from  the  6000°K  energy  distribution  to  the  7000°K 
energy  distribution.  For  burst  heights  greater  than  25  km, 
linear  interpolation  is  used  from  a  7000°K  distribution  at 
25  km  to  an  8000°K  distribution  at  35  km.  Above  35  km,  the 
8000°K  distribution  is  used.  To  increase  the  flexibility  of 
LOWTHRM ,  this  routine  was  adopted  as  a  default  option.  The 
user  may  either  specify  the  blackbody  temperature  under  con¬ 
sideration  or  defer  to  the  computed  default  value. 

The  thermal  yield  of  each  device  under  consideration  is 
determined  through  use  of  an  empirical  two-dimensional  linear 
interpolation  routine.  The  routine  was  developed  by  the 
Analysis  Division,  Air  Force  Weapons  Laboratory  (AFWL) , 
Kirtland  AFB,  New  Mexico.  It  is  based  on  the  latest  Sputter 
data  with  thermal  efficiency  calculated  as  a  function  of 
height  of  burst  and  thermal  yield  as  a  function  of  total 
weapon  yield.  The  thermal  efficiency,  T,  is  given  by 

T  =  exp  ((-.358  -  .009H  +  7.14  x  10-4HZ)  x  2.30)  x 
exp  ((-1.25  x  10_5H3  +  6.42  x  10-8H4)  x  2.30) 

where  H  is  the  height  of  burst  in  kilometers.  The  thermal 
yield  is  then  expressed  by  the  relation 

Y  =  (T  x  W)  x  1012  (cal)  (18) 

where  Y  is  the  thermal  yield,  T  is  the  thermal  efficiency, 
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and  W  is  the  total  weapon  yield  in  kilotons.  Because  of  the 
atmospheric  model  limitations,  the  altitudes  used  in  the 
routine  should  be  limited  to  100  kilometers.  The  yield  input 
has  a  recommended  upper  limit  of  100  megatons  (Ref.  12:28). 

Direct  Fluence  Model 

Models  by  Johnson  and  others  had  used  discrete  energy 
bands  to  divide  the  spectrum  under  consideration.  From  these 
bands,  average  transmittance  values  were  determined  and  then 
used  across  the  entire  band.  From  transmittance  studies  made 
while  examining  the  basic  LOWTRAN5  code,  it  was  found  that 
this  broad  band  approximation  was  unnecessary.  The  basic 
code  provided  transmissivity  data  from  350  cm-1  to  40000  cm-1 
(0.25  -  28.5  microns)  at  20  cm-1  frequency  increments.  This 
fine  band  structure  was  retained  and  improves  the  accuracy 
of  the  code. 

Computed  transmissivity  at  each  frequency  interval  is 
weighted  with  a  Planckian  spectrum  based  on  the  source  tem¬ 
perature  which  is  input  by  the  user  or  computed  by  Johnson's 
interpolation  routine.  The  LOWTHRM  code  begins  calculations 
at  350  cm”1  (28.5  microns)  and  increments  in  20  cm-1  steps 
to  40000  cm-1  (0.25  microns).  To  initialize  the  first  weight¬ 
ing  factor  to  28.5  microns,  a  routine  was  incorporated  to 
generate  a  Planckian  integral  to  this  frequency  before  trans¬ 
mittance  calculations  were  begun.  This  integral  was  approxi¬ 
mated  by  a  trapezoidal  sum  of  the  form: 


22 


B (350) 


(19) 


15 


350 

Z 

u=5 


u3 (Au) 
eu  -  1 


where  B  is  the  value  of  the  integral  sum  at  28.5  microns, 
u  =  h  =  Planck's  constant  (6.6756  x  10“27erg  -  sec), 
v  =  frequency  (sec  -1),  K  =  Boltzmann's  constant  (1.38054  x 
10” 1 6erg/°K) ,  T  =  source  temperature  (°K) ,  and  Au  =  As 

the  fluence  is  computed  for  each  frequency  interval  throughout 
the  range  of  the  program,  this  integral  is  continually  updated 
so  that  the  weight  for  each  increment  is  correct  for  the 
source  temperature  in  use.  The  weight  for  each  band  is  then: 

BW  (v)  =  B (v)  -  B(v  -  Av)  (20) 


Adding  conversion  factors  to  convert  kilometers  to  centi¬ 
meters  and  kilotons  to  calories  produces  output  fluences  in 
units  of  calories  per  square  centimeter.  The  direct  fluence 
at  frequency  v  can  then  be  approximated  by: 


Q(v) 


( BW ( V ) )  (Y)  (T(V))  (1012) 

(4ir)  (r  x  10s)  2 


(cal/cm2) 


(21) 


where  Q  is  the  direct  fluence  at  frequency  v,  BW  is  the  band 
weight,  Y  is  the  thermal  yield,  T  is  the  transmittance,  and 
r  is  the  path  length  (km) .  The  sum  over  all  wavelengths  is 
the  total  direct  fluence  on  target  and  is  approximated  by 

40000 

Q.  =  Z  Q ( v)  (cal/cm2)  (22) 

*■  v=350 

where  v  is  incremented  in  20  cm“J  steps. 
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Scattered  Fluence  Model 


Research  revealed  two  possible  models  to  approximate 
scattered  fluence.  The  Air  Force  Weapons  Laboratory  (AFWL) 

TRAP  Code  (Ref.  11)  uses  an  empirical  buildup  factor  to 
account  for  diffuse  scattering  by  the  atmosphere.  To  use 
this  buildup  factor,  scattering  coefficients  based  on  height 
and  wavelength  had  to  be  computed  for  each  frequency  interval. 
LOWTHRM  was  modified  to  use  this  buildup  factor  and  results 
were  compared  with  published  values.  The  computed  scattered 
fluences  were  much  lower  than  expected  and  the  buildup  factor 
model  was  not  retained  in  LOWTHRM. 

Tedesco  (Ref.  24)  developed  a  single-scatter  integral 
for  use  as  a  computer  model  in  thermal  fluence  calculations. 

This  model  assumes  that  any  radiation  which  is  scattered  more 
than  once  along  its  path  will  not  be  incident  upon  the  receiver 
in  the  time  interval  under  consideration.  Due  to  the  long 
mean  free  paths  of  thermal  radiation  in  air,  this  assumption 
is  a  reasonable  approximation.  A  further  assumption  is  that 
the  transmission  medium  is  undisturbed  by  the  nuclear  radiation 
flux  from  the  detonation.  A  close  examination  of  his  work  re¬ 
vealed  that  his  basic  transmission  model  closely  parallelled 
the  1962  U.S.  Standard  Atmosphere  and  Rural  Aerosol  Model 
options  available  in  LOWTHRM.  A  telephone  conversation  with 
Mr.  Frank  Kneizys  of  the  Air  Force  Geophysics  Laboratory  (AFGL) , 
Hanscom  AFB,  Massachusetts  verified  that  all  necessary  values 
for  scattering  coefficients  and  absorber  amounts  were  available 
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through  the  L0WTRAN5  code.  The  single-scatter  integral  was 
incorporated  into  LOWTHRM.  Because  the  model  now  required 
double  integrations  over  all  scattering  planes  and  frequencies, 
execution  time  was  a  major  problem.  Tedesco  had  complained 
of  this  even  over  the  fourteen  bands  that  his  model  used. 

The  complication  was  much  worse  in  LOWTHRM  since  there  are 
approximately  two  thousand  increments  in  the  range  of  the 
program.  Results  close  to  those  published  by  Tedesco  could 
be  obtained,  but  run  times  were  well  in  excess  of  four  hun¬ 
dred  seconds  execution  time.  The  range  of  LOWTHRM  was  reduced 
to  0.25  to  4.0  microns;  making  the  assumption  that  very  little 
of  the  spectrum  was  beyond  this  wavelength.  This  change  did 
produce  minor  reductions  in  the  direct  fluences.  However, 
the  execution  time  was  still  in  the  100  to  300  second  range 
depending  upon  the  particular  problem  and  number  of  scattering 
levels  involved. 

This  execution  time  problem  is  resolved  in  LOWTHRM  by 
modifications  to  the  scatter  integral.  Tedesco  had  used  a 
fifth  degree,  least  square  polynomial  curve  fit  to  approxi¬ 
mate  Mie  and  Rayleigh  scattering.  Curve  fit  coefficients 
varied  as  the  frequency  increased  throughout  the  range  of 
his  program  (Ref.  24:84).  This  curve  fit  was  discarded  and 
is  replaced  in  LOWTHRM  by  simple  isotropic  scatter.  The 
single-scatter  integral  was  rewritten  in  terms  of  angular 
integrations  only  and  then  normalized  to  the  direct  fluence. 

By  normalizing  to  the  direct  fluence,  the  scatter  integral 
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produces  a  ratio  which  is  independent  of  frequency.  By  mul¬ 
tiplying  this  ratio  times  the  direct  fluence  at  any  wavelength, 
the  scattered  amount  for  that  wavelength  is  then  obtained. 

Thus,  the  requirement  for  scatter  integrations  at  each  fre¬ 
quency  increment  is  eliminated  and  greatly  reduces  execution 
time.  The  ratio  is  developed  in  the  following  manner. 

The  scatter  geometry  is  shown  in  Figure  2.  Distances 
have  been  normalized  to  the  vertical  displacement  of  the 
burst  point  from  the  scattering  plane.  The  outward  unit  sur¬ 
face  normal  of  the  receiver  has  components  Xn,  Yn,  and  ZR. 

Xn  is  in  the  horizontal  plane  parallel  to  the  scattering  plane 
and  is  positive  toward  the  burst  point.  is  in  the  hori¬ 

zontal  plane  and  normal  to  Xn»  and,  due  to  symmetry,  may  be 
considered  positive  in  any  direction.  Zn  is  normal  to  the 
scattering  plane  and  is  positive  toward  the  scattering  plane. 

The  elemental  area  dA  (normalized  to  burst  height  squared) 
is 

dA  =  tan0sec20d0d<j>  (23) 

The  radiation  normally  incident  at  the  elemental  area  dA  is 
T2<3c^s  30/4tt.  Parameter  q  is  the  total  radiant  energy  inten¬ 
sity  from  the  source  and  T2  is  the  transmittance  from  the 
burst  to  dA.  The  fraction  of  this  radiation  scattered  to  the 
receiver  point  is  T-^cos  2  3/4nH2  where  T^  is  the  transmittance 
from  dA  to  the  receiver.  The  direct  radiation  at  the  receiver 
point  is  qT^/ 4 tt  ( (H— 1)  2  +  G2).  The  ratio  of  the  scattered 
radiation  to  the  direct  is  then: 
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Fig.  2.  Scatter  Geometry  (Ref.  11:49) 


If  jWJ.B  MlIfKII 


scattered 
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scattered 
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4 it  (  (H-l)  2+G2)  T-qcos 3 6T,cos 2 0tan8 sec2 0d0d<{> 

- 4  4“ - ^ - <24> 


qT, 


(4tt )  (4ttH2) 
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T^4ir 


'4  4 
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H‘ 


scattered  T0T, ( (H-l) 2+G2 )  cos2  0sin0d8d<{> 

-  =  - /  f  -  (26) 

direct  T^tt  ^  H2 


The  component  of  scattered  radiation  normal  to  the  receiver 
must  now  be  determined.  The  unit  surface  normal  may  be  writ¬ 
ten  in  vector  form  as  X  I  +  Y  ]  +  Z  ic.  A  unit  vector  from 

n  n  n 

the  receiver  to  the  area  dA  is: 

(HtanBcosyi  +  Htangsinyj  +  Hk)cos0/H 
The  dot  product  of  these  two  vectors  is  the  cosine  of  the 
angle  between  the  surface  normal  and  the  ray  from  the  area 
dA.  Multiplying  the  integrand  of  Equation  (26)  by  this  dot 
product  then  yields  the  component  of  the  scattered  radiation 
normal  to  the  receiver.  So,  the  ratio  of  the  scattered  to  the 
direct  may  be  written  as: 


4  4 


T2T- ( (H-l) Z+G2) 

RATIO  =  -=—Z -  X 

T^4tt 

cos 3  $sin0  (Xntan0cosy+Yntan£5siny+Zn)  d0d<f> 
”  ~  ~  H2  ~  '  ~ 


(27) 


Before  the  integration  can  be  performed,  0  and  y  must  be  ex- 

l 

pressed  in  terms  of  8  and  <J>.  Referring  again  to  Figure  2. 

HtanBcosy  +  tan0cos4>  =  G  (28) 
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tangcosy  = 


G  -  tan0cos<{> 


H 


Htangsiny  =  tan0sin<{> 
tan0sin<J) 


tanfSsiny  = 


H 


H 


cos  g  = 


(H2+G2+tan2  0-2Gtan0cos<(>)  ^ 

Substituting  the  above  equations  into  Equation  5, 

T-T- ( (H-l) 2+G2)  H3sin0 

RATIO  =  -=— = - L  L 


®  ^  (H2+G2+tan20-2Gtan0cos<f))  ^ 
Collecting  terms  we  have. 


T^4ir 


fQ  f<t> 


T-T- ( (H-l) Z+GZ)H 

RATIO  =  -  X 

T^47r 

sin0  (X  (G-tan0cos<^)+Y  (tan0sin<J>)  +Z  H)d0d<f> 


n 


n 


n 


(H2+G2+tanz0-2Gtan0cos({)) 
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(29) 

(30) 

(31) 

(32) 


(33) 


(34) 


All  values  are  known  except  for  T2  and  T^ •  These  are  approx¬ 
imated  by  the  expression: 


Z  1  N 

/  T(Z)dZ=±  E  T(Z)AZ 

0  Ni=l 


(35) 


where  AZ  is  taken  as  1  kilometer  for  ease  of  calculation  and 
N  is  the  distance  from  the  source  to  the  scattering  plane  for 
T2  and  the  distance  from  scattering  plane  to  the  receiver 
for  Tj.  Values  for  transmittance  at  altitude  Z,  T(Z),  were 
determined  from  the  LOWTRAN5  code  for  each  model  atmosphere 
by  taking  average  transmittances  for  each  one  kilometer  layer 
of  the  atmosphere  up  to  a  height  of  one  hundred  kilometers. 
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This  was  accomplished  for  both  the  23  km  and  the  5  km  visi¬ 
bility  profiles.  Transmittances  to  two  hundred  kilometers 
were  approximated  as  equal  to  the  one  hundred  kilometer  value. 
These  values  were  stored  in  data  arrays  for  use  in  subsequent 
scattering  calculations  when  called  by  the  sort  routine. 

The  sort  routine  determines  the  number  of  scattering 
planes  for  each  particular  problem.  Three  special  cases  are 
considered  and  depicted  in  Figure  3.  In  case  one  the  number 
of  scattering  planes  is  limited  to  twice  the  number  of  kilo¬ 
meters  in  the  path  length,  i.e.  20  planes  for  a  10  km  path. 

In  this  case  the  burst  height  is  exactly  equal  to  the  path 
length.  Case  two  considers  path  lengths  greater  than  burst 
altitude.  In  this  case  the  number  of  scattering  planes  is 
limited  to  the  number  of  kilometers  in  the  path  length  plus 
the  burst  height,  i.e.  17  for  a  path  length  of  10  km  and  burst 
height  of  7  km.  The  final  case  is  for  burst  height  greater 
than  path  length.  Scattering  plane  determination  is  the  same 
as  for  case  two,  i.e.  17  for  a  10  km  burst  and  7  km  path 
length.  The  use  of  the  path  length  to  limit  the  number  of 
scattering  planes  is  based  on  studies  of  the  scattering  ratio. 
Contributions  to  scatter  rapidly  decrease  as  the  distance  from 
burst  to  scatter  plane  increases.  The  scattering  ratio  for 
planes  at  distances  greater  than  the  path  length  produced 
values  that  were  about  one  per  cent  of  the  contributions  from 
planes  close  to  the  burst.  Because  these  contributions  are 
so  small,  the  sort  routine  limits  the  number  of  planes  as  des¬ 
cribed  in  this  paragraph. 
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To  actually  perform  the  required  angular  integrations, 
a  subroutine  used  in  the  AFWL  Trap  Code  (Ref.  11)  was  modified 
for  use  in  LOWTHRM.  Integrations  over  each  scattering  plane 
are  performed  as  called  by  the  sort  routine.  As  the  direct 
fluence  at  each  wavelength  is  determined  and  summed,  the 
scattered  fluence  is  calculated  and  added  to  its  sum.  The 
approximation  is: 

40000 

Q  =  (RATIO)  E  Q  (v)  (36) 

S  v=350 

where  Qg  is  the  total  scattered  fluence,  RATIO  is  the  com¬ 
puted  ratio  of  scattered  to  direct  fluence,  Q(v)  is  the 
direct  fluence  at  frequency  v,  and  v  is  incremented  in  20  cm-1 
steps. 

Ground  Reflection  Model 

In  order  to  be  a  complete  fluence  model,  LOWTHRM  includes 
a  ground  reflected  contribution  with  variable  surface  albedo. 
The  albedo  values  along  with  corresponding  reflecting  sur¬ 
faces  are  listed  in  Table  IV  (Ref.  4).  The  ground  reflected 
fluence  was  treated  similarly  to  the  scattered  fluence  in 
that  a  normalized  reflection  integral  was  developed.  The 
problem  is  simplified  somewhat  since  there  is  only  one  plane 
required  for  ground  reflection  instead  of  the  several  requi-ed 
for  scatter  computations.  One  additional  assumption  is  that 
the  ground  reflectance  can  be  approximated  by  Lambert's  cosine 
law.  This  results  in  an  additional  factor  of  2cos6  in  the 
numerator  of  Equation  (24).  The  integral  is  developed  in 
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the  same  manner  as  the  scatter  integral  and  will  not  be  re¬ 
peated  here.  Adding  the  variable  albedo  constant,  the  final 
reflection  integral  is: 

T.Tc-eHC  (H-l)  2+G2) 

RR  =  - -  x 

T^tt 

(37) 

sin0  (Xn  (G-tan0cos<f)  +Yntan0sin$+ZnH)  d0d<J> 

fe  ^ 

(H2+G2+tan2  8-2Gtan8cos<{>)  2 

where  RR  is  the  reflection  ratio  of  the  reflected  to  the 
direct  fluence,  is  the  transmittance  from  the  ground  to 
the  receiver,  T^  is  the  transmittance  from  the  burst  to  the 
receiver,  and  e  is  the  albedo  constant.  Angles  0  and  <{>  and 
variables  G  and  H  are  specified  in  Figure  4.  As  in  the  scat¬ 
ter  approximation,  all  values  of  the  reflection  integral  are 
known  except  T4  and  Tj.  These  are  approximated  by  the 
expression: 

Z  1  N 

f~T  (Z)  dZ“  Z  T(Z)AZ  (38) 

Ni=l 

where  AZ  is  1  km,  T(Z)  is  the  average  transmittance  at  alti¬ 
tude  Z,  and  N  is  the  distance  from  the  source  to  the  ground 
for  and  the  distance  from  the  ground  to  the  receiver  for 


With  the  completion  of  the  ground  reflected  component, 
the  model  has  been  completely  specified.  The  only  possible 
addition  would  be  to  include  a  variable  cloud  albedo  and 
scattering  approximations  within  cloud  layers.  Since  no 
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Table  IV 

Variable  Surface  Albedos 


SURFACE 

ALBEDO 

Desert 

0.310 

Various  Fields 

0.190 

Green  Forest 

0.115 

Grass 

0.305 

Bare  Ground 

0.185 

Dry  Sand 

0.230 

Wet  Sand 

0.140 

Snow  or  Ice 

0.680 

Rough  Water 

0.205 

Average  Water 

0.100 

Shock-Frothed  Water 


0.600 


V.  Results  and  Conclusions 

As  mentioned  in  the  introduction,  the  purpose  of  this 
thesis  is  to  incorporate  atmospheric  transmissivity  factors 
into  a  new  thermal  radiation  code.  To  determine  the  validity 
of  this  model,  results  are  compared  against  Tedesco's  Single- 
Scatter  model,  (Ref. 24)  the  AFWL  Snare  code  (Ref. 11),  Johnson's 
Virgin  Fluence  model  (Ref. 12),  and  the  Thermal  Routine  from 
Horizon  Technology,  Inc.  which  has  been  developed  for  use  on 
the  Texas  Instruments  TI-59  calculator  (Ref.  8) .  Table  V 
shows  a  comparison  for  LOWTHRM  against  Tedesco's  Single- 
Scatter  model  and  the  AFWL  Snare  code.  All  fluences  are  for 
horizontal  paths  at  the  altitudes  indicated.  The  LOWTHERM 
model  atmosphere  used  was  the  1962  U.S.  Standard  Atmosphere 
with  the  Rural  (23km)  Aerosol  Model.  Since  Tedesco  used  a 
ground  reflection  albedo  of  0.35,  LOWTHRM  was  modified  to 
accept  that  value  for  comparative  purposes.  A  comparison  of 
the  values  listed  in  Table  V  shows  strong  agreement  up  to  an 
altitude  of  twenty-five  kilometers.  For  the  runs  at  twenty- 
five  and  thirty-five  kilometers,  LOWTHRM  yields  higher  flu¬ 
ences  on  target.  This  may  be  due  to  the  better  transmission 
calculations  LOWTHRM  offers,  a  difference  in  the  spectrum  used 
to  weight  each  frequency  interval,  or  a  combination  of  both. 

The  results  for  two  available  runs  of  the  AFWL  Snare  code  are 
listed  for  comparative  purposes  of  current  Air  Force  esti¬ 
mates  of  thermal  fluence.  It  should  be  noted  that  these 
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estimates  are  based  upon  a  visibility  of  ten  miles  at  sea 
level,  while  LOWTHRM  is  based  upon  the  more  accurate  parameter 
of  meteorological  range. 

Table  V 

Fluence  Comparison  for  LOWTHRM,  Single-Scatter 
Model,  and  AFWL  Snare  Codes 


HEIGHT 

(km) 

PATH 

(km) 

YIELD 

(kt) 

LOWTHRM 
(cal/cm2 ) 

SINGLE- SCATTER 
(cal/cm2 ) 

SNARE 
(cal/cm2 ) 

2 

10 

1000 

25.17 

25.49 

28.5 

5 

10 

1000 

35.41 

35.77 

- 

7 

10 

1000 

35.99 

36.18 

36.5 

15 

8 

100 

5.34 

5.41 

- 

25 

8 

100 

5.69 

4.79 

- 

35 

16 

1000 

15.85 

13.09 

- 

Table  VI  is  a  comparison  between  Johnson's  Virgin  Fluence 
model,  LOWTHRM,  and  the  Thermal  Routine  from  the  Nuclear  Wea- 
ns  Effects  Programs  Weapons  Manual. 


Table  VI 

Fluence  Comparison  for  Virgin  Fluence 
Model,  LOWTHRM,  and  HTI  Thermal  Routine 


TARGET 

BURST 

RANGE 

YIELD 

VIRGIN 

LOWTHERM 

THERMAL 

(km) 

(km) 

(km) 

(kt) 

(cal/cm2 ) 

(cal/cm2 ) 

(cal/cm2 ) 

2 

5 

200 

13.21 

18.32 

19.26 

5 

10 

200 

3.04 

6.03 

4.53 

10 

12 

1000 

12.20 

25.82 

17.23 

i> 


It  should  be  noted  that  Johnson's  Virgin  Fluence  Model  is  for 
direct  fluence  values  only.  The  values  listed  in  Table  VI 
differed  from  LOWTHRM  direct  fluences  only  by  seven  to  ten 
percent.  The  Thermal  Routine  by  Horizon  Technology,  Inc. 
appears  to  be  a  "data  fit"  of  some  sort  which  yields  fluences 
fairly  close  to  LOWTHRM  and  the  other  models  investigated  in 
this  report. 

In  conclusion,  LOWTHRM  is  a  fast  running  code  that  can 
calculate  thermal  fluence  for  a  variety  of  model  atmospheres 
and  aerosol  conditions.  The  "exact"  atmospheric  transmissivity 
factors  used  in  the  direct  fluence  calculations  provide  the 
basis  for  a  high  degree  of  accuracy  while  the  adoption  of  the 
normalized  scatter  integral  greatly  reduces  execution  times 
with  no  loss  in  output  fluences. 
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Appendix  A 
LOWTHRM  User 1 s  Guide 

LOWTHRM  is  activated  by  submission  of  one  data  card. 

The  data  card  is  read  in  a  free  field  format  of  the  form 

READ*, MODEL, IHAZE, ITYPE, ISEASN, VIS,W1A,T1A, IPOSIT, IABEDO, 

HI, H2, RANGE 

Parameter  one,  MODEL,  is  used  to  specify  one  of  the  following 

atmospheres  by  entering  an  integer  value  between  0  and  7. 

0  =  Horizontal  path  and  meteorological  data  are  user 
specified 

1  =  Tropical  Atmosphere 

2  =  Midlatitude  Summer 

3  =  Midlatitude  Winter 

4  =  Subarctic  Summer 

5  =  Subarctic  Winter 

6  =  1962  U.S.  Standard  Atmosphere 

7  =  New  Model  is  user  specified  using  radiosonde  data 

IHAZE  is  an  integer  value  entered  to  select  the  appropriate 
aerosol  mode. 

1  =  Rural  (23km) 

2  =  Rural  (5km) 

3  =  Maritime  (23km) 

4  =  Maritime  (5km) 

5  =  Urban  (5km) 
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6  =  Tropospheric  (50km) 

7  =  User  Defined 

8  =  Fogl  (0.2km) 

9  -  Fog2  (0.5km) 

Parameter  three,  ITYPE,  is  an  integer  used  to  indicate  the 
type  of  atmospheric  path. 

1  =  Horizontal  path 

2  =  Vertical  or  slant  path  between  two  altitudes 

ISEASN  specifies  seasonal  visibility  profiles. 

1  =  Spring-Summer 

2  =  Fall-Winter 

Parameter  five  is  VIS.  This  is  a  real  value  which  specifies 
the  visual  range  (km)  at  sea  level.  Any  value  other  than  zero 
will  activate  the  interpolation/extrapolation  routine  con¬ 
tained  in  LOWTHRM.  If  VIS  is  set  equal  to  zero,  a  default 
value  is  assigned  equal  to  the  visibility  specified  by  I HAZE. 

W1A  is  a  real  value  which  is  the  weapon  yield  in  kilotons. 

T1A  is  the  temperature  (°K)  of  the  Planckian  spectrum  used  to 
weight  the  transmission  calculations.  If  the  user  sets  T1A 
equal  to  zero,  LOWTHRM  computes  a  default  value  based  on  the 
height  of  burst.  TlA  is  entered  as  a  real  value. 

Parameter  eight  is  IPOSIT.  It  is  used  to  specify  the  re¬ 
ceiver  orientation  and  is  entered  as  a  single  integer  value. 
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IPOSIT  specifies  the  direction  cosines  of  the  unit  vector 
normal  to  the  target  surface.  Values  of  X  are  positive  in 
the  direction  of  the  burst.  Values  of  Z  are  positive  towards 
the  ground.  If  direction  cosines  other  than  those  listed 
below  are  required,  the  user  can  select  his  own  values  by 
setting  IPOSIT  equal  to  0  and  then  entering  the  desired  values, 
i.e. - , TlA, 0 , . 866, .866, . 866 , IALBEDO , Hi ,H2 , RANGE. 


X 

y 

Z 

1 

= 

1. 

0. 

0. 

2 

= 

.866 

0. 

-.5 

3 

= 

.5 

0. 

-.866 

4 

= 

-.5 

0. 

-.866 

5 

= 

-.866 

0. 

-.5 

6 

=- 

■1. 

0. 

0. 

7 

= 

-.866 

0. 

.5 

8 

= 

-.5 

0. 

.866 

9 

= 

0. 

0. 

1. 

10 

= 

.5 

0. 

.866 

11 

= 

.866 

0. 

.5 

12 

= 

.866 

.5 

0. 

13 

= 

.5 

.866 

0. 

14 

= 

0. 

1. 

0. 

15 

= 

-.5 

.866 

0. 

16 

= 

-.866 

.5 

0. 

17 

= 

0. 

.5 

-.866 

18 

= 

0. 

.  866 

-.5 

19 

= 

0. 

.866 

.5 

20 

as 

0. 

.5 

.866 

Parameter  nine  is  IALBEDO.  This  is  an  integer  value  used  to 
select  the  ground  reflection  constant.  The  user  can  select 
his  own  ground  reflection  constant  by  setting  IALBEDO  to  0 

43 


followed  by  the  desired  constant,  i.e. - IPOSIT,0, . 725, Hi, - . 

1  =  Desert  (0.310) 

2  =  Various  Fields  (0.190) 

3  =  Green  Forest  (.115) 

4  =  Grass  (0.305) 

5  =  Bare  Ground  (0.185) 

6  =  Dry  Sand  (0.280) 

7  =  Wet  Sand  (0.140) 

8  =  Snow  or  Ice  (0.680) 

9  =  Rough  Water  (0.205) 

10  =  Average  Water  (0.100) 

11  =  Shock-Frothed  Water  (0.600) 

Parameter  ten  is  Hi.  This  is  a  real  value  which  corresponds 
to  the  target  altitude  in  kilometers. 

H2  is  the  burst  attitude  in  kilometers. 

The  final  parameter  is  RANGE.  This  specifies  the  slant  range 
between  the  burst  and  receiver.  It  is  a  real  value  with  units 
of  kilometers. 

To  facilitate  use  of  LOWTHRM,  it  is  recommended  that  the 

binary  source  deck  be  stored  on  a  permanent  file,  LWTHERM, 

cycle  one.  The  data  deck  should  be  stored  on  permanent  file, 

LWTHERM,  cycle  two.  The  program  may  then  be  attached  and 

run  as  follows: 

xxx, T30 , CM100000 , STANY. 

ATTACH , WES , LWTHERM , CY= 1 
ATTACH , TAPE  5 , LWTHERM , CY=  2 
WES. 

7/8/9 
DATA  CARD 
7/8/9 
6/7/8/9 
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The  above  procedure  eliminates  the  need  to  read  the  source 
program  and  data  deck  for  each  run. 

There  are  two  special  cases  which  may  be  of  interest  to 
some  users.  Case  one  is  for  a  horizontal  path  when  meteorolo¬ 
gical  data  is  used.  In  this  instance,  MODEL  is  set  to  zero. 
After  the  data  card,  enter  a  second  card  with  the  following 
information: 

Hi  (km) ,  pressure  (mb) ,  temperature  (degrees  C) ,  dew 
point  temperature  (deg.  C) ,  percent  relative  humidity,  water 
vapor  density  (gm/m3),  ozone  density  (gm/m3),  RANGE  (km). 

When  using  this  option.  Format  statement  429  in  Subroutine 

NSMDL  specifies  the  input  to  be  used  for  the  above. 

Case  two  allows  the  user  to  build  a  new  model  atmosphere. 
On  the  data  card  set  MODEL  to  seven  and  add  data  cards  for 
each  altitude  level  in  the  new  model  atmosphere.  Each  new 
card  should  contain  the  following: 

altitude  (km) ,  pressure  (mb) ,  temperature  (Deg.  C) ,  Dew 
point  (Deg.  C) ,  relative  humidity,  water  vapor  density  (gm/m3), 
ozone  density  (gm/m3),  aerosol  number  density  (cm-3),  VIS, 

I HAZE,  ISEASN,  and  IVULCN. 


Format  435  from  Subroutine  NSMDL  is  used  for  case  two. 

Other  parameters  within  LOWTHRM  which  may  be  varied  by 
the  user  include  Ml,  M2,  M3,  RO,  TBOUND,  IVULCN,  ANGLE,  BETA, 
IM,  ML,  LEN,  and  Y.  With  the  exception  of  Y,  these  para¬ 
meters  are  normally  set  equal  to  zero  for  operation  of  the 
program,  but  may  be  varied  by  the  user.  To  change  them,  a 
card  substitution  is  required  in  the  source  program  for  each 
parameter  that  is  to  be  changed.  This  means  the  user  must 
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remove  the  card  that  states  M1=0  and  replace  it  with  the  option 
chosen,  i.e.  Ml=l.  Each  parameter  is  annotated  with  comment 
cards  in  the  listing  so  that  each  is  easy  to  find  and  substi¬ 
tution  is  simple  when  required. 

Ml  is  an  integer  used  to  vary  temperature  altitude  profiles. 

0  =  Normal  Operation 

1  =  Tropical 

2  =  Midlatitude  Summer 

3  =  Midlatitude  Winter 

4  =  Subarctic  Summer 

5  =  Subarctic  Winter 

6  =  1962  U.S.  Standard 

M2  is  an  integer  used  to  vary  water  vapor  altitude  profiles. 

0  =  Normal  operation 

1  =  Tropical 

2  =  Midlatitude  Summer 

3  =  Midlatitude  Winter 

4  =  Subarctic  Summer 

5  =  Subarctic  Winter 

6  =  1962  U.S.  Standard 

M3  is  used  to  vary  ozone  altitude  profiles. 

0  =  Normal  Operation 

1  =  Tropical 

2  =  Midlatitude  Summer 


3  =  Midlatitude  Winter 
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4  =  Subarctic  Summer 


5  =  Subarctic  Winter 

6  =  1962  U.S.  Standard 

RO  is  a  real  value  which  specifies  the  radius  of  the  earth 
at  the  particular  location  at  which  the  calculation  is  to  be 
performed.  If  RO  is  set  equal  to  zero,  default  value  is  the 
earth  radius  for  the  model  atmosphere  selected. 

TBOUND  is  the  earth  temperature  in  degrees  Kelvin.  TBOUND 
equals  zero  defaults  to  the  air  temperature  of  the  model 
atmosphere  selected. 

IVULCN  is  an  integer  defining  the  ten  to  thirty  kilometer 
aerosol  visibility  profile. 

0  =  Default  to  Stratospheric  background 

1  =  Stratospheric  background 

2  =  Aged  Volcanic  type/Moderate  content 

3  =  Fresh  Volcanic  type/High  content 

4  =  Aged  Volcanic  type/High  content 

5  =  Fresh  Volcanic  type/Moderate  content 

ANGLE  is  the  zenith  angle  at  HI.  If  ANGLE  is  set  equal  to 
zero,  a  default  value  is  computed  from  values  specified  for 
HI,  H2,  and  RANGE. 

BETA  is  the  earth  center  angle  s>xbtended  by  Hi  and  H2.  If 
BETA  is  set  equal  to  zero,  a  default  value  is  computed  from 
values  specified  for  HI,  H2,  and  RANGE. 
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IM  is  an  integer  value  for  use  with  MODEL  equal  to  seven. 
This  corresponds  to  special  case  two  previously  mentioned 
and  IM  should  be  set  equal  to  one. 


ML  is  an  integer  value  for  the  number  of  atmospheric  levels 
to  be  read  in  for  special  case  two.  For  horizontal  paths 
ML=1.  For  slant  paths  ML  equals  the  number  of  one  kilometer 
levels  traversed. 

LEN  is  a  path  length  control  parameter.  For  some  calculations 
in  which  ANGLE  or  BETA  have  been  specified,  LOWTHRM  may 
determine  that  there  are  two  possible  paths  for  computations. 
In  this  event,  a  message  will  indicate  the  shorter  path  was 
computed.  To  select  the  longer  path,  set  LEN  equal  to  one. 

The  last  parameter  to  be  discussed  is  Y,  the  thermal  yield. 
Some  literature  specifies  methods  for  computing  the  thermal 
yield  other  than  the  model  used  in  LOWTHRM.  To  change  the 
thermal  yield  approximation,  the  user  may  substitute  source 
program  cards  at  the  point  documented  in  the  program  listing. 
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Appendix  B 
Sample  Problems 

This  appendix  contains  six  sample  problems  to  further 
aid  the  user  in  running  LOWTHRM.  These  problems  are  not  "all 
inclusive"  in  exercising  the  options  made  available.  They  do, 
however,  give  the  potential  user  output  for  specific  problems 
and  an  insight  into  what  can  be  done  with  LOWTHRM. 
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Problem  1 

The  first  sample  problem  is  to  compute  the  ^otal  fluence 
on  target  for  a  200  kiloton  weapon  with  the  receiver  at  1 
kilometer  height  and  the  burst  at  3  kilometer  height.  Slant 
range  is  4  kilometers.  The  1962  U.S.  Standard  Atmosphere  is 
used  with  the  Rural  Aerosol  model.  Sea  level  visibility  is 
23  kilometers.  Ground  reflection  is  from  desert  and  receiver 
direction  cosines  are  X  =  1.0,  Y  =  0.,  and  Z  =  0.  The  season 
is  spring. 

The  data  card  parameters  for  problem  one  are: 

MODEL  =  6 
I HAZE  =  1 
ITYPE  =  2 
I SEAS N  =  1 

VIS  =  0.  (Defaults  to  I HAZE  value) 

WlA  =  200. 

TlA  =  0.  (Defaults  to  a  value  based  on  height  of  burst) 

I POSIT  =  1 
I ALBEDO  =  1 
HI  =  1. 

H2  =  3. 

RANGE  =  4. 

Since  the  data  card  is  read  in  a  free  field  format,  it  may  be 
punched  as: 

6, 1,2, 1,0. ,200. ,0. ,1,1,1. ,3. ,4. 

Output  is  as  follows: 


50 


INPUT  PARAMETERS . . . 

WEAPON  YIELD  IS  200.  KILOTONS 

SLANT  PATH  BETWEEN  ALTITUDES  Hi  and  H2  WHERE  Hi  =  1. 
KM  H2  =  3.000  KM 

PATH  LENGTH  IS  4.00  KM 

HAZE  MODEL  =  23.0  KM  VISUAL  RANGE  AT  SEA  LEVEL 

MODEL  ATMOSPHERE  6  =  1962  US  STANDARD 

HAZE  MODEL  1  =  RURAL  VIS  =  23.0KM 

SEASON  =  SPRIG  SUMM 

GROUND  REFLECTION  IS  FROM  DESERT 

RECEIVER  DIRECTION  COSINES  ARE  X  =  1.000  y  =  0.000 
Z  =  0.000 

VERTICAL  PROFILE  AEROSOL  MODEL  =  STRAT  BKGR 

OUTPUT  PARAMETERS . . . 

AT  THE  SOURCE . . . 

IR  THERMAL  YIELD  IS  40.53  KILOTONS 
VISIBLE  THERMAL  YIELD  IS  31.45  KILOTONS 
UV  THERMAL  YIELD  IS  10.65  KILOTONS 
TOTAL  THERMAL  YIELD  IS  83.70  KILOTONS 

AT  THE  RECEIVER. . . 

IR  UNREACTED  FLUENCE  IS  14.88  CAL/ CM2 
VISIBLE  UNREACTED  FLUENCE  IS  12.60  CAL/ CM 2 
UV  UNREACTED  FLUENCF  IS  2.79  CAL/CM2 
TOTAL  UNREACTED  FLUENCE  IS  30.27  CAL/CM2 
SCATTERED  FLUENCE  IS  5.73  CAL/CM2 
REFLECTED  FLUENCE  IS  2.55  CAL/CM2 
TOTAL  FLUENCE  ON  TARGET  IS  38.55  CAL/CM2 
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Problem  2 

Problem  number  two  is  to  repeat  problem  one  but  change 
the  receiver  direction  cosines  to  X  =  0,  Y  =  0,  and  Z  =  1. 

The  only  parameter  different  than  problem  number  one  is 
IPOS IT. 

IPOS IT  =  9 

The  Data  card  is  now: 

6, 1, 2, 1,0., 200., 0., 9, 1,1., 3. ,4. 

The  output  on  the  following  page  reflects  the  change  in 
receiver  orientation: 
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INPUT  PARAMETERS . . 

WEAPON  YIELD  IS  200.  KILOTONS 

SLANT  PATH  BETWEEN  ALTITUDES  HI  AND  H2  WHERE  HI  =  1 
KM  H2  =  3.000  KM 

PATH  LENGTH  IS  4.000  KM 

HAZE  MODEL  =  23.0  KM  VISUAL  RANGE  AT  SEA  LEVEL 

MODEL  ATMOSPHERE  6  =  1962  US  STANDARD 

HAZE  MODEL  1  =  RURAL  VIS  =  23.0KM 

SEASON  =  SPRIG  SUMM 

GROUND  REFLECTION  IS  FROM  DESERT 

RECEIVER  DIRECTION  COSINES  ARE  X=  0.000  Y  =  0.000 

z  =  1.000 

VERTICAL  PROFILE  AEROSOL  MODEL  =  STRAT  BKGR 

OUTPUT  PARAMETERS . . . 

AT  THE  SOURCE . . . 

IR  THERMAL  YIELD  IS  40.53  KILOTONS 
VISIBLE  THERMAL  YIELD  IS  31.45  KILOTONS 
UV  THERMAL  YIELD  IS  10.65  KILOTONS 
TOTAL  THERMAL  YIELD  IS  83.70  KILOTONS 
AT  THE  RECEIVER... 

IR  UNREACTED  FLUENCE  IS  0.00  CAL/ CM2 
VISIBLE  UNREACTED  FLUENCE  IS  0.00  CAL/ CM2 
UV  UNREACTED  FLUENCE  IS  0.00  CAL/ CM2 
TOTAL  UNREACTED  FLUENCE  IS  0.00  CAL/CM2 
SCATTERED  FLUENCE  IS  7.56  CAL/ CM2 
REFLECTED  FLUENCE  IS  3.61  CAL/CM2 
TOTAL  FLUENCE  ON  TARGET  IS  11.16  CAL/ CM2 
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Problem  3 


The  third  example  repeats  problem  one  with  the  receiver 
at  20  kilometers,  the  burst  at  30  kilometers,  and  a  slant 
range  of  10  kilometers. 

HI  =  20. 

H2  =  30. 

RANGE  =  10. 

All  other  parameters  are  the  same  as  for  problem  one. 

The  Data  card  is  now: 

6, 1,2, 1,0. ,200. ,0, 1,1, 20. ,30. , 10. 

The  output  on  the  following  page  indicates  an  increase 
in  thermal  yield  due  to  a  change  in  the  height  of  burst  and 
a  decrease  in  fluence  due  to  increased  path  length. 
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INPUT  PARAMETERS . . . 

WEAPON  YIELD  IS  200.  KILOTONS 

SLANT  PATH  BETWEEN  ALTITUDES  HI  AND  H2  WHERE  HI  =  20 
KM  H2  =  30.000  KM 

PATH  LENGTH  IS  10.000  KM 

HAZE  MODEL  =  23.0  KM  VISUAL  RANGE  AT  SEA  LEVEL 

MODEL  ATMOSPHERE  6  =  1962  US  STANDARD 

HAZE  MODEL  1  =  RURAL  VIS  =  23.0KM 

SEASON  =  SPRIG  SUMM 

GROUND  REFLECTION  IS  FROM  DESERT 

RECEIVER  DIRECTION  COSINES  ARE  X  =  1.000  Y  =  0.000 

Z  =  0.000 

VERTICAL  PROFILE  AEROSOL  MODEL  =  STRAT  BKGR 
OUTPUT  PARAMETERS . . . 

AT  THE  SOURCE . . . 

IR  THERMAL  YIELD  IS  36.43  KILOTONS 
VISIBLE  THERMAL  YIELD  IS  42.28  KILOTONS 
UV  THERMAL  YIELD  IS  24.28  KILOTONS 
TOTAL  THERMAL  YIELD  IS  108.26  KILOTONS 
AT  THE  RECEIVER. . . 

IR  UNREACTED  FLUENCE  IS  2.88  CAL/CM2 
VISIBLE  UNREACTED  FLUENCE  IS  3.31  CAL/ CM2 
UV  UNREACTED  FLUENCE  IS  1.27  CAL/CM2 
TOTAL  UNREACTED  FLUENCE  IS  7.46  CAL/CM2 
SCATTERED  FLUENCE  IS  .16  CAL/CM2 
REFLECTED  FLUENCE  IS  .03  CAL/CM2 
TOTAL  FLUENCE  ON  TARGET  IS  7.65  CAL/ CM2 


Problem  4 

* 

•v  The  fourth  example  problem  determines  the  thermal  fluence 

on  target  for  a  1  megaton  weapon  with  the  target  and  burst 
co-altitude  at  7  kilometers.  The  path  length  is  10  kilometers. 
The  model  is  Tropical  with  the  Rural  Aerosol  Model.  Sea  level 
visibility  is  23  kilometers.  The  season  is  spring  with  ground 
reflection  user  defined  as  0.35.  Receiver  direction  cosines 
are  user  defined  as  X  =  0.866,  Y  =  0.866,  Z  =  0.866. 

The  data  card  contains  the  following: 

MODEL  =  1 
I HAZE  =  1 
ITYPE  =  1 
ISEASN  =  1 
VIS  =  0. 

4  ■ 

WlA  «  1000. 

TlA  =  0. 

IPOSIT  =  0, . 866 , . 866 , . 866  (0  specifies  user  defined 

followed  by  input  values) 

I ALBEDO  =  0,.35  (0  specifies  user  defined 

followed  by  input  value) 

HI  =  7. 

H2  =  7. 

RANGE  =  10. 

Output  is  on  the  following  page. 


\  ■ 
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INPUT  PARAMETERS . . . 

WEAPON  YIELD  IS  1000.  KILOTONS 

HORIZONTAL  PATH,  ALTITUDE  =  7.000  KM,  PATH  LENGTH  = 
10.000  KM 

HAZE  MODEL  =  23.0  KM  VISUAL  RANGE  AT  SEA  LEVEL 
MODEL  ATMOSPHERE  1  =. TROPICAL 
HAZE  MODEL  1  =  RURAL  VIS  =  23.0  KM 
SEASON  =  SPRIG  SUMM 

GROUND  REFLECTION  IS  FROM  USER  DEFIN 

RECEIVER  DIRECTION  COSINES  ARE  X  =  0.866  Y  =  0.866 
Z  =  0.866 

VERTICAL  PROFILE  AEROSOL  MODEL  =  STRAT  BKGR 

OUTPUT  PARAMETERS... 

AT  THE  SOURCE . . . 

IR  THERMAL  YIELD  IS  197.81  KILOTONS 
VISIBLE  THERMAL  YIELD  IS  153.49  KILOTONS 
UV  THERMAL  YIELD  IS  51.98  KILOTONS 
TOTAL  THERMAL  YIELD  IS  408.53  KILOTONS 
AT  THE  RECEIVER. . . 

IR  UNREACTED  FLUENCE  IS  11.33  CAL/CM2 
VISIBLE  UNREACTED  FLUENCE  IS  9.22  CAL/CM2 
UV  UNREACTED  FLUENCE  IS  1.91  CAL/ CM2 
TOTAL  UNREACTED  FLUENCE  IS  22.46  CAL/ CM2 
SCATTERED  FLUENCE  IS  7.84  CAL/ CM2 
REFLECTED  FLUENCE  IS  3.83  CAL/CM2 
TOTAL  FLUENCE  ON  TARGET  IS  34.12  CAL/CM2 
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Problem  5 

The  fifth  problem  is  to  compute  the  thermal  fluence  on 
target  for  a  1  megaton  weapon  with  the  receiver  on  the  ground 
and  the  burst  altitute  at  4  kilometers.  The  slant  range  is 
8  kilometers.  The  1962  U.S.  Standard  Atmosphere  is  to  be 
used  with  the  Urban  Aerosol  Model.  Visual  range  at  sea  level 
is  5  kilometers.  The  season  is  summer  and  ground  reflection 
is  from  bare  ground.  Receiver  direction  cosines  are  X  =  .5, 

Y  =  0.,  and  Z  =  .866. 

Data  card  parameters  are: 

MODEL  =  6 

I HAZE  =  5 

1TYPE  =  2 

ISEASN  =  1 

VIS  =  0. 

W1A  =  1000. 

TlA  =  0 . 

I POSIT  =  10 

I ALBEDO  =  5 

HI  =  0. 

H2  =  4. 

RANGE  =  8 . 

Output  is  on  the  following  page. 
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INPUT  PARAMETERS . . . 

WEAPON  YIELD  IS  1000.  KILOTONS 

SLANT  PATH  BETWEEN  ALTITUDES  Hi  AND  H2  WHERE  HI  =  0.000 
KM  H2  =  4.00  KM 

PATH  LENGTH  IS  8.000  KM 

HAZE  MODEL  =  5 . 0  KM  VISUAL  RANGE  AT  SEA  LEVEL 
MODEL  ATMOSPHERE  6  =  1962  US  STANDARD 
HAZE  MODEL  5  =  URBAN  VIS  =  5.0  KM 
SEASON  =  SPRIG  SUMM 

GROUND  REFLECTION  IS  FROM  BARE  GRND 

RECEIVER  DIRECTION  COSINES  ARE  X  =  .500  Y  =  0.000 
Z  =  .866 

VERTICAL  PROFILE  AEROSOL  MODEL  =  STRAT  BKGR 
OUTPUT  PARAMETERS . . . 

AT  THE  SOURCE . . . 

IR  THERMAL  YIELD  IS  200.66  KILOTONS 
VISIBLE  THERMAL  YIELD  IS  155.71  KILOTONS 
UV  THERMAL  YIELD  IS  52.73  KILOTONS 
TOTAL  THERMAL  YIELD  IS  414.42  KILOTONS 
AT  THE  RECEIVER. . . 

IR  UNREACTED  FLUENCE  IS  3.33  CAL/ CM 2 
VISIBLE  UNREACTED  FLUENCE  IS  1.03  CAL/ CM2 
UV  UNREACTED  FLUENCE  IS  .07  CAL/ CM2 
TOTAL  UNREACTED  FLUENCE  IS  4.43  CAL/ CM2 
SCATTERED  FLUENCE  IS  2.61  CAL/CM2 
REFLECTED  FLUENCE  IS  0.00  CAL/CM2 
TOTAL  FLUENCE  ON  TARGET  IS  7.04  CAL/CM2 
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Problem  6 


...  rA 


The  sixth  and  last  example  problem  illustrates  special 
case  one  mentioned  in  the  LOWTHRM  User's  Guide.  In  this 
problem,  specific  meteorological  data  is  used  for  a  horizontal 
path  at  10  kilometer  height.  The  weapon  yield  is  1  megaton 
with  a  path  length  of  10  kilometers.  The  Rural  Aerosol  Model 
is  used  with  a  sea  level  visibility  of  50  kilometers.  The 
season  is  fall.  The  receiver  direction  cosines  are  X  =  1., 

Y  =  0,  and  Z  =  0.  Ground  reflection  is  from  desert.  Pressure 
is  1000  millibars,  ambient  temperature  is  10°C,  relative 
humidity  is  40%,  and  the  Midlatitude  Winter  Ozone  profile  is 
to  be  used. 

The  data  card  parameters  are: 

MODEL  =  0 

I HAZE  =  1 

ITYPE  =  1 

ISEASN  =  2 

VIS  =  50. 

W1A  =  1000. 

T1A  =  0. 

IPOS IT  =  1 

I ALBEDO  =  1 

HI  =  10. 

H2  =  10. 

RANGE  =  10. 
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A  second  data  card  is  used  for  this  problem.  It  contains 
values  for  target  height  (km) ,  pressure  (mb) ,  temperature 
(Deg.  C) ,  Dew  Point  temperature  (Deg.  C)  ,  percent  relative 
humidity,  water  vapor  density  (gm/m3),  ozone  density  (gm/m3) , 
and  path  length  (km) .  In  this  case  the  parameters  are  as 
follows  using  format  429  found  in  Subroutine  NSMDL: 

Target  height  =  10. 

Pressure  =  1000. 

Temperature  =  10. 

Dew  Point  temperature  =  0.  (Default  value  is  computed) 

Percent  relative  humidity  =  40. 

Water  Vapor  Density  =  0 .  (Not  used  in  this  problem) 

Ozone  Density  =  0.  (Midlatitude  Winter  profile  used — 

see  below) 

Path  length  =  10. 

Since  the  Midlatitude  Winter  Ozone  profile  is  to  be  used, 
a  source  program  card  substitution  is  required.  The  user  must 
remove  the  card  "M3  =  0"  from  the  program  and  add  "M3  =  3"  to 
load  the  appropriate  values. 

Output  is  essentially  the  same  as  before.  The  line 
specifying  model  atmosphere  selected  is  not  printed  when 
MODEL  is  set  to  zero.  For  this  problem  the  output  is  on  the 
following  page. 
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INPUT  PARAMETERS . . . 

WEAPON  YIELD  IS  1000.  KILOTONS 

HORIZONTAL  PATH,  ALTITUDE  =  10.000  KM,  PATH  LENGTH  = 

10.000  KM 

HAZE  MODEL  =  50.0  KM  VISUAL  RANGE  AT  SEA  LEVEL 

HAZE  MODEL  1  =  RURAL  VIS  =  50.0KM 

SEASON  =  FALL  WINTR 
GROUND  REFLECTION  IS  FROM  DESERT 

RECEIVER  DIRECTION  COSINES  ARE  X  =  1.000  Y  =  0.000 

Z  =  0.000 

VERTICAL  PROFILE  AEROSOL  MODEL  =  STRAT  BKGR 

OUTPUT  PARAMETERS . . . 

AT  THE  SOURCE . . . 

IR  THERMAL  YIELD  IS  198.82  KILOTONS 
VISIBLE  THERMAL  YIELD  IS  154.28  KILOTONS 
UV  THERMAL  YIELD  IS  52.24  KILOTONS 
TOTAL  THERMAL  YIELD  410.63 
AT  THE  RECEIVER. . . 

IR  UNREACTED  FLUENCE  IS  11.29  CAL/CM2 
VISIBLE  UNREACTED  FLUENCE  IS  10.45  CAL/ CM2 
UV  UNREACTED  FLUENCE  IS  1.61  CAL/CM2 
TOTAL  UNREACTED  FLUENCE  IS  23.35  CAL/ CM2 
SCATTERED  FLUENCE  IS  7.49  CAL/CM2 
REFLECTED  FLUENCE  IS  .83  CAL/CM2 
TOTAL  FLUENCE  ON  TARGET  IS  31.67  CAL/ CM2 
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\  Fortran  computer  program  LOWTHRM  is  described  for  calculating 
nuclear  thermal  fluence  incident  upon  a  target  area.  Atmospheric 
transmissivity  factors  in  the  spectral  region  0.25  to  28.5  microns 
are  determined  through  use  of  the  LOWTRAN5  computer  code.  The 
program  provides  a  choice  of  six  model  atmospheres  covering  seasonal 
and  latitudinal  variations  from  sea  level  to  100  km,  eight  haze 
models,  and  accounts  for  molecular  absorption,  molecular  scattering, 
and  aerosol  extinction.  Atmospheric  refraction,  earth  curvature 
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effects,  thermal  scattering,  and  thermal  ground  reflection 
contributions  are  included. 
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